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The von Neumann architecture for a classical computer comprises a central processing unit and a memory 
holding instructions and data. We demonstrate a quantum central processing unit that exchanges data with 
a quantum random-access memory integrated on a chip, with instructions stored on a classical computer. 
We test our quantum machine by executing codes that involve seven quantum elements: Two superconduct- 
ing qubits coupled through a quantum bus, two quantum memories, and two zeroing registers. Two vital 
algorithms for quantum computing are demonstrated, the quantum Fourier transform, with 66 % process 
fidelity, and the three-qubit Toffoli OR phase gate, with 98 % phase fidelity. Our results, in combination 
especially with longer qubit coherence, illustrate a potentially viable approach to factoring numbers and 
implementing simple quantum error correction codes. 

Quantum processorsG^ based on nuclear magnetic resonance!^, trapped ions^ '", and semiconducting de- 
viceslJJ^were used to realize Shor's quantum factoring algorithmic and quantum error correction!^. The quantum 
operations underlying these algorithms include two-qubit gates^E the quantum Fourier transformO^l, and three- 
qubit Toffoli gatesli^l^. In addition to a quantum processor, a second critical element for a quantum machine is a 
quantum memory, which has been demonstrated, e.g., using optical systems to map photonic entanglement into 
and out of atomic ensembles'*^. 

Superconducting quantum circuitslSlhave met a number of milestones, including demonstrations of two-qubit 
g^)-gg|5ii5j |i7 | i9 | 20 | ^jjjj jjjg advanced control of both qubit and photonic quantum statesl^ESESlll demonstrate 
a superconducting integrated circuit that combines a processor, executing the quantum Fourier transform and a 
three-qubit Toffoli-class OR gate, with a memory and a zeroing register in a single device. This combination of a 
quantum central processing unit (quCPU) and a quantum random-access memory (quRAM), which comprise two 
key elements of a classical von Neumann architecture, defines our quantum von Neumann architecture. 

In our architecture (Fig. lA), the quCPU performs one-, two-, and three-qubit gates that process quantum 
information, and the adjacent quRAM allows quantum information to be written, read out, and zeroed. The 
quCPU includes two superconducting phase qubits'^^L^'^ Q ^ and Q2, connected through a coupling bus provided 
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Fi gure 1: The quantum von Neumann architecture. (A) The quCPU (blue box) includes two qubits Q-j^ and Q2 and the bus resonator B. 
The quRAM (magenta boxes) comprises two memories Mj^ and and two zeroing registers Zj^ and Zj. The horizontal dotted lines indicate 
connections between computational elements. The vertical direction represents frequency, where the memory and zeroing registers are fixed 
in frequency, while the qubit transition frequencies can be tuned via z-pulses (grey dashed double aiTows). (B) Swap spectroscopy' for Qj^ 
(left) and Qj (right): Qubit excited state |e) probability (color scale) vs. z-pulse amplitude (vertical axis) and delay time At (horizontal 
axis), after exciting the qubit with a 7r-pulse. At zero z-pulse amplitude the qubits are at their idle points, where they have an energy relaxation 



time T , ~ 400 ns. A separate Ramsey experiment yields the qubits' dephasing time T, 



deph 



200 ns. By tuning the z-pulse amplitude, the 



qubit transition frequencies /g and /q can be varied between ~ 5.5 and 8 GHz. For z-pulse amplitudes indicated by B and Mj^ for , and 

by B and for Qj, the "chevron pattern" of a qubit-resonator interaction is observed'. The transition frequencies of B, Mj^, and Mj are 
/g = 6.82 GHz, = 6.29 GHz, and = 6.34 GHz, respectively. From the chevron oscillation we obtain the qubit-resonator coupling 

strengths, which for both the resonator bus and the memories are ~ 20MHz (splitting) for the |g) -M- |e) qubit transition, and \/2 faster for 
the |e) |f) transition (|g), |e), and |f) are the three lowest qubit states)^^. For all resonators T^^^ ~ 4 /xs. Swap spectroscopy also reveals 
that the qubits interact with several modes associated with spurious two-level systems. Two of them, Z-^ and Zj, are used as zeroing registers. 
Their transition frequencies are = 6.08 GHz and =7.51 GHz, respectively, with coupling strength to the qubits of ~ 17 MHz. 



by a superconducting microwave resonator B. The quRAM comprises two superconducting resonators Mj and 
M2 that serve as quantum memories, as well as a pair of zeroing registers Zi and Zj, two-level systems that are 
used to dump quantum information. The chip geometry is similar to that in Refs.121221 with the addition of the two 
zeroing registers. Figure IB shows the characterization of the device by means of swap spectroscopy . 

The computational capability of our architecture is displayed in Fig. 2A, where a 7-channel quantum circuit, 
yielding a 128 dimensional Hilbeit space, executes a prototypical algorithm. First, we create a Bell state between 
Qi and Q2 using a seiies of vr-pulse, -\/iSWAP, and iSWAP operations (step I, a to c)^^. The corresponding density 
matrix p^j [Fig. 2C (I)] is measured by quantum state tomography. The Bell state is then written into the quantum 
memories and M2 by an iSWAP pulse (step 11)123 leaving the qubits in their ground state |g), with density 
matrix [Fig. 2C (II)]. While storing the first Bell state in and M2, a second Bell state with density matrix 
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Figure 2: Programming the quantum von Neumann architecture. (A) Quantum algorithm comprising 7 independent channels interacting 
tiu'ough five computational steps. Dotted and solid lines represent channels in the ground and excited/superposition states, respectively. A 
black rectangle represents a 7r-pulse; two crosses connected by a solid line a \/ iSWAP; an open and a closed circle connected by a single 
arrow an iSWAP; oblique arrows indicate decay from a zeroing register. (B) Calibration of the zeroing gates. Each qubit is prepared in |e), 
interacts on resonance with its zeroing register for a time r^, and its probability measured, with plotted vs. (large and small blue 
circles). The solid green line is a decaying cosine fit to the data. The black arrows indicate the zeroing time for each qubit. (C) Density 
matrices Pyj, Pju), • ■ ■ , P^y) of Qi"Q2 ^'^'^ ^™ ench step in A (scale key on bottom left). Grey arrows: Ideal state. Red and black arrows 
and black dots: Measured state (black arrows indicate errors). The off-diagonal elements of p^j, p^jny and p^yj have different angles because 



of dynamic phases^". Fidelities: 



(I) 



■^(V) ~ 0.606 it 0.003. Concurrences: 



0.772 ± 0.003, JPy,) = 
0.593 ± 0.006, C„ 



0.916 ± 0.002, JP'yjjj = 
- 0.029 ± 0.005, C(ii,) : 



0.913 ± 0.002, and 
= 0.019 ± 0.005, and 



0.689 ± 0.003, JP(,vj = 
- ^ v.v.., - 0.436 ± 0.007, C(,v) ^ 

CjYj = 0.345 lb 0.008. (D) Comparison of fidelity as a function of storage time r^, for a Bell state stored in Qj and (blue circles) vs. 
that stored in Mj^ and (magenta squares; error bars smaller than symbols). The solid lines are exponential fits to data. (E) As in D, but for 
the concurrence C. In D and E the vertical black dotted line indicates the time delay (~ 59 ns) associated with memory storage, with respect 
to storage in the qubits, due to the writing and reading operations (II) and (V) in A. 



/'(Ill) F^S- 2C (m)] is created between the qubits, using a sequence similar to the first operation (step III, a to c). 

In order to re-use the qubits and Qj, for example to read out the quantum information stored in the mem- 
ories Ml and M2, the second Bell state has to be dumped^S. This is accomplished using two zeroing gates, by 
bringing on resonance with Zj^ and Q2 with for a zeroing time r^, corresponding to a full iSWAP (step IV). 
Figure 2B shows the corresponding dynamics, where each qubit, initially in the excited state |e), is measured in 
the ground state |g) after ~ 30 ns. The density matrix p^y) of the zeroed two-qubit system is shown in Fig. 2C 
(IV). Once zeroed, the qubits can be used to read the memories (step V), allowing us to verify that, at the end 
of the algorithm, the stored state is still entangled. This is clearly demonstrated by the density matrix shown in 
Fig. 2C (V). 

The ability to store entanglement in the memories, which are characterized by much longer coherence times 
than the qubits, is key to the quantum von Neumann architecture. We demonstrate this capability in Fig. 2, D and 
E, where the fidelity and concurrence metricsl^of the Bell states stored in Mj^ and M2 are compared to those for 
the same states stored in Qj^ and Q2. The experiment is performed as in Fig. 2 A, but eliminating steps (III) and 
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Figure 3: The quantum Fourier transform. (A) (Left) Quantum logic circuit of a CZ-4> gate (enclosed in a grey box) for IQ1Q2) = |ee). 
The |f) state of Qj is indicated by a dashed line. The process where Qj^ acquires the phase is represented by a pair of open/closed circles, 
connected by a single arrow in an arc shape. All other symbols are as in Fig. 2A. (Right) Shorthand symbol for the CZ-(j> gate. Although the 
gate unitary matrix is symmetric, the symbol shows the asymmetric implementation of the gate. (B) Time-domain swaps between the states 
iQjB) = |el) and |fO), where we plot the probability (color scale) vs. interaction time At and detuning 5q g. The solid black Hne 

indicates combinations of interaction time and detuning that completely depopulate the non-computational |f) state. The three black dots on 
this line correspond to a CZ-tt, CZ-7r/2, and CZ-0.28 gate (see far right). The fourth black dot (outside the line) corresponds to a i/2CZ-7r 
gate (see bottom-left), where the |e) state has been shelved to the non-computational |f) state. (C) Phase (f> acquired by Qj^ as a function 
of (5q g. The blue dots indicate experimental data and the solid green line the theory of Eq. flP*. (D) Fidelity (blue symbols) and 
EOF (magenta "x" symbols) of measured density matrices vs. 5^ g. (E) (Left to Right) EJensity matrices = jg, '^t' 
obtained when <f> = 0.28, <j> = tt/2, and (p = Trrad in Eq.[2](scale key on bottom left). The arrows are color-coded as in Fig. 2C. The 
measured fidelities are js = O-^Sl ± 0.064, T^^^ = "-735 ± 0.017, and T„ = 0.741 ± 0.030, and EOF are S^ ^s = 0-020 ± 0.055 
(lower bound Sq 2a — 0)" ~ 0-^06 it 0.031, and = 0.401 it 0.062. (F) (Top-Left) Logic circuit for a two-qubit quantum Fourier 

transform and, (Bottom), real part of the corresponding Xm matrixHS. The process fidelity for the real and imaginary (not shown) part of Xm 
is J-"^ = 0.657 it 0.014. The confidence intervals are estimated from 10 measurements for pg jg, 6 for p^^2 ^^'^ Pn' ^^d 15 for Xm- 

(IV). For the qubits, the storage time r^, is defined as the wait time at the end of step (I), prior to measuring the 
qubit states, whereas for the resonators the wait time is that between the write and read steps. The fidelity of the 
qubit states decays to below 0.2 after 400 ns, while for the states stored in the memories it remains above 0.4 up to 
~ 1.5 /is. Most importantly, after only 100 ns the state stored in the qubits does not preserve any entanglement, as 
indicated by a zero concurrence, whereas the memories retain their entanglement for at least 1.5 fis (Fig. 2E). We 
expect taking advantage of our architecture in long computations, where qubit states can be protected and reused 
by writing them into, and reading them out of, the long-lived quRAM. 

Two-qubit universal gates are a vital resource for the operation of the quCPU*^. A variety of such gates 
have been implemented in superconducting circuitsl^ISEIIllIol vvith some recent demonstrations of quantum al- 
gorithms^E^. Control Z-tt (CZ-tt) gates are readily realizable with superconducting qubits, due to easy access to 
the third energy state of the qubit, effectively operating the qubit as a qutrit)3^^El2l] However, CZ-tt gates are just 
a subset of the more general class of CZ-4> gates, obtained for the special case where the phase = tt. In our 
architecture, the full class of CZ-(/) gates, with 4> from ~ to tt, can be generated by coupling a qutrit close to 
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resonance with a bus resonator. 

Figure 3A shows the quantum logic circuit that generates the CZ-(/) gate (Left) and a shorthand symbol for 
the gate (Right). The logic circuit demonstrates the nontrivial case where qubits and Q2 are brought from 
their initial ground state to IQ1Q2) = |ee) by applying a tt -pulse to each qubit. The excitation in Q2 is then 
transferred into bus resonator B, and Qi's |e) ^ |f) transition brought close to resonance with B for the time 
required for a 27r-rotation, where the states iQj^B) = |el) and |fO) are detuned by a frequency 5q g, which we 
term a "semi-resonant condition." In this process Qi acquires the phasel^^ 



where g is the coupling frequency between |el) and |fO). The final step is to move the excitation from B back 
into Q2. 

The time-domain swaps of IQ^B) between the states |el) and |fO) are shown in Fig. 3B, where the solid black 
line indicates the detunings and corresponding interaction times used to generate any phase < ^ tt (ideally 
0^-0 when (5q g 00). These phases are measured by performing two Ramsey experiments on for each 
value of the detuning Sq^^, one with B in the |0) state, and the other with B in the |1) state. The relative phase 
between the Ramsey fringes corresponds to the value of (j) for the CZ-(j) gatel^, as shown in Fig. 3C. 

A more sophisticated version of this experiment is performed by initializing Qj^ and Q2 each in the superposi- 
tion state |g) + |e). We move Qj's state into B, perform a CZ-0 gate with < (j> ^ tt, move the state in B back into 
Q2, rotate Qj^'s resulting state by tt/2 about the y-axis, and perform a joint measurement of Qi and Q2. Ideally, 
this protocol permits to create two-qubit states ranging from a product state for = to a maximally-entangled 
state for 4> = tt. In the two-qubit basis set — {|gg)i l^g), |ge), |ee)}, the general density matrix of such 
two-qubit states reads 



/O \ 

1/2 {l-e~"t')/4: {l + e-"t')/4: 

{l-e"t>)/4: (l-cos0)/4 (-isin0)/4 

\0 (l + e*'^)/4 (zsin</.)/4 (l + cos0)/4y 



(2) 



Figure 3D shows the fidelity and entanglement of formation (EOF)Elof two-qubit states generated using 70 values 
of (f). Figure 3E shows three examples of for (f) = 0.28, = tt/2, and (p = tt, respectively. 

The state generated using (p = tt/2 plays a central role in the implementation of the two-qubit quantum 
Fourier transform. Neglecting bit-order reversal, the quantum Fourier transform can be realized by applying a 
Hadamard gate to Q2, followed by a CZ-7r/2 gate between Qj^ and Q2, and finally a Hadamard on Qj^l2l212[ as 
sketched in Fig. 3F (Top-Left). Representing the input state of the transform as \x) (position) and the output as 
\p) (momentum), assuming \x) E M2 ^i^^ the indexes x and p are integers, with p E {0,1,2, 3}, the output state 
\p) — J2x=o '^^ ^^^^ corresponding to a 4 x 4 unitary operator This operator can be fully characterized by 

means of quantum process tomography's^] which allows us to obtain the Xm matrix"^ shown in Fig. 3F (Bottom). 

Finally, by combining the CZ-0 and zeroing gates, we can implement a Toffoli-class gateE^EEII the three- 
qubit OR phase gate. This gate, combined with single qubit rotations, is sufficient for universal computation. 
A Toffoli gate is a doubly-controlled quantum operation, where a unitary operation is applied to a target qubit 
subject to the state of two control qubits. The canonical Toffoli is a doubly-controlled NOT gate; here we consider 
a doubly-controlled phase gate, which is equivalent through a change of basis of the target qubit. In the canonical 
Toffoli gate, the control gate is applied if both control qubits, AND Q2, are in state |e). In our case, the 
control gate is applied conditionally if the controls Qi OR Q2 are in |e). Additionally, we have implemented a 
three-qubit gate for the logical function XOR, which, even though not a Toffoli-class gate, helps to understand the 
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Figure 4: Three-qubit gates: The XOR phase gate and the Toffoli-class M gate. (A) Quantum logic circuit for the XOR pliase gate. (B) 
(Top) XOR-gate truth table. (Bottom) Ramsey fringes associated with the truth table, showing the probability of measuring Qj in |e), vs. 
the Ramsey phase ip, for the control input states in A^2- Black and magenta dots: phase. Blue and green dots: vr phase. The solid lines are 
least-squares fits to the data used to extract the truth-table phases. (C) Quantum phase tomography for the XOR gate: Phase (p^i^^-^ — </'|5go) ■ 
for each state \lmn) £ M^. Black open boxes: Ideal values. Pink areas: Measured values with corresponding confidence intervals (black 
lines). (D) Quantum logic circuit for the M gate, implemented as a 1/2 CZ-tt gate (cf. Fig. 3B) between Qj^ and B (half-dot/half-open circle 
connected by solid line), followed by a CZ-tt gate between Qj and B, and a second 1/2 CZ-tt gate between Qj^ and B. The dotted black lines 
connecting the two 1/2 CZ-tt gates indicate qubit shelving to the |f) state. (E) As in panel B, but for the M gate. (F) As in C, but for the M 
gate. 



more complex OR gate. 

The quantum logic circuits for the XOR and OR gates are drawn in Fig. 4, A and D. The control qubits are Qi 
and Q2 and the target is the bus resonator B, effectively acting as the third qubit (as only the states |0) and |1) of 
B are used). The XOR gate is realized as a series of two CZ-tt gates between the controls and the target, and the 
OR gate as the series 1/2 CZ-tt, CZ-tt, and 1/2 CZ-tt, in an "M-shape" configuration. 

The truth table for the XOR gate is displayed in Fig. 4B (Top). The control qubits and Q2 are assumed 
to be in one of the states in A^2' while the target B is in |0) + The target acquires a phase tt, corresponding 
to a "true" result, only when the controls are in the state IQ1Q2) = |ge) or |eg). For the other non-trivial case 
IQ1Q2) — k^)' the target acquires phase, corresponding to a "false" result. This is due to the action of the two 
CZ-TT gates, giving a global phase tt when either of the controls is in |e), and a phase 2tt (equivalent to a phase) 
when both are in |e). 

The truth table can be experimentally measured by performing Ramsey experiments on the target, one for each 
pair of control states. The experiments are realized by, (i), preparing Q2 in the superposition state |g) + |e) by 
means of a TT/2-pulse; (ii), moving the state from Q2 into B, thus creating a |0) + |1) state in B; (Hi), preparing Qi 
and Q2 in each possible pair of control states in A^2 by means of TT-pulses; (iv), performing the XOR gate; (v), 
zeroing Q2 into Z2 at the end of the XOR gate; (vi), moving the final target state from B into the zeroed Q2; (vii), 
completing the Ramsey sequence on Q2 with a second TT/2-pulse with variable rotation axis relative to the pulse 
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in ( i). The measurement outcomes are displayed in Fig. 4B (Bottom), together with the least-squares fits used to 
extract the phase information associated with each value of the truth table. The Ramsey fringes for the two control 
states |ge) and |eg) are inverted relative to the reference state |gg), as expected from the XOR gate truth table. 

In general, given the Q1-Q2-B basis set A^g = {|ggO), |ggl), |geO), |gel), |egO), 
|egl), |eeO), |eel)}, the vector r™"^ of the diagonal elements associated with the ideal unitary matrix of the XOR 
gate reads 

r^""^ = (1, 1, 1, -1, 1, -1, 1, 1) , (3) 

while all off-diagonal elements of the matrix are zero. Each element t^°^ can be expressed as a complex expo- 
nential e* "^i'™"> , with \lmn) e A4^. The phase can be either 0, when t^^^ — 1, or tt, when t^°^ = — 1. 
Among the eight values of 4)^1,^^^^ only seven are physically independent, as the element e*"*ii!so> can be factored, 
reducing the set of possible phases to 0|;^„^ — </'|ggo)' "^^^^ \lmn) e Ai^ — {|ggO)}. 

In analogy to the truth-table for the target B, a table with four phase differences can also be obtained for 
the controls and Q2, resulting in a total of twelve phase differences. These differences can be measured by 
performing Ramsey experiments both on the target and the control qubits. It can be shown that from the twelve 
phase differences, one can obtain the seven independent phases associated with the diagonal elements r^'-"*'^, thus 
realizing a quantum phase tomography of the Toffoli gatel^Sl Figure 4C displays the phase tomography results for 
our experimental implementation of the XOR gate. 

The truth table associated with the M gate is reported in Fig. 4E (Top), where the only difference from the XOR 
gate is the phase tt acquired by the target B when the controls Qi and Q2 are loaded in state IQ1Q2) = |ee). In this 
case, the action of the first 1/2 CZ-tt gate between Qj^ and B shelves the |1) state from B to the non-computational 
state |f) in Q^, where it remains until the second 1/2 CZ-tt gate. Moving the state of Qj^ outside the computational 
space during the intermediate CZ-tt gate between Q2 and B effectively turns off the CZ-tt gate '^'^^'. The target 
B thus only acquires a total phase tt due to the combined action of the two 1/2 CZ-tt gates (cf. Fig. 4D). The 
experimental truth table obtained from Ramsey fringes is shown in Fig. 4E (Bottom). 

The vector of the diagonal elements associated with the ideal unitary matrix of the M gate is = 
(1, 1, 1, —1, 1, —1, 1, —1). A similar procedure as for the XOR gate allows us to obtain the quantum phase 
tomography of the M gate (Fig. 4F). 

Quantum phase tomography makes it possible to define the phase fidelity of the XOR and M gate, 

-F^^l-^, (4) 

where is the gate root-mean-square phase error, with an upper bound of tt. For the XOR gate we find that 
= 0.954 ± 0.004, and for the M gate = 0.979 ± 0.003. 

Our results provide optimism for the near-term implementation of a larger-scale quantum processop^*^ based 
on superconducing circuits. Our architecture shows that proof-of-concept factorization algorithmsl^^E and simple 
quantum error correction 

codesiasMI 

might be achievable using this approach. 
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Materials and Methods 
Statistical errors 

In this section, we analyze the statistical properties of the experimental data shown in the main text. First, we 
explain how to simulate statistical errors. This procedure was used to estimate the confidence intervals for the data 
of Fig. 2 in the main text. Second, we describe how statistical errors were obtained from statistical ensembles of 
independent measurements. This procedure was used for the data of Fig. 3 in the main text. Third, we discuss the 
estimation of statistical errors due to fits to the data. This procedure was used for the data of Fig. 4 in the main 
text. 

Simulation of statistical errors 

In this subsection, we discuss two important sources of statistical errors in our data: Errors associated with qubit's 
measurement (binomial-type errors) and errors due to jitter/fluctuations in the electronics (phase errors). Assuming 
binomial-type and phase errors, we describe the procedures used to simulate the confidence intervals for the 
elements and metrics of the density matrices shown in Fig. 2C of the main text. 

( i) Binomial-type errors are inherent to our qubit measurement process, where the measurement is repeated a 
fixed number of times N, each measurement trial has two possible outcomes, i.e., qubit being in the ground 
state |g) with probability p„ or in the excited state |e) with probability — 1 — p^, the probability p^ is to 
good approximation the same for each trial, and the trials can be considered to be statistically independent. 
The measurement outcome associated with |g) is counted as 0, and that associated with |e) as 1. Under 
these assumptions, the qubit measurement process can be described by a binomial distribution. 



12 



Given a statistical sample consisting of N measurement outcomes (i.e., a statistical sample from 
a Bernoulli distribution with parameter p^), the maximum likelihood estimator of (i.e., the estimated 
probability) is given by 

fc=l 

where X^ represents the fc-th outcome among the N measured. There are several ways to compute a 
confidence interval for the parameter p^. The most common result is based on the approximation of the 
binomial distribution with a normal distribution. This represents a good approximation in our experiments, 
where the number of measurements N is large (typically N ^ 1500). In this case, it can be shown that a 
confidence interval for the parameter p^ is given by 



^'e ± 2(l_a/2) Y ' ^ ' =P,± ^(l-«/2) : (S6) 

where ^(i_q,/2) is the (1 — a/2) percentile of a standard normal distribution. For example, for a 0.95 (95%) 
confidence interval, we set a = 0.05, so that ~ 1.96. When analyzing our data we approximate 

the percentile 1.96 with 2, thus obtaining a slightly wider confidence interval; 

( ii) Phase errors are mostly due to the phase jitter/fluctuations in the room-temperature cables and electronics 
used to measure the qubits. In order to quantify such errors, the following experiment was performed. First, 
we initialized one of the two qubits, e.g., qubit Q^, in the ground state, \Qi) — |g); second, we applied to 
Qi a 7r/2 unitary rotation about the y-axis, Ry^^, bringing the qubit into the state \Qi) — (|g) + |e))/ \/2. 
This state is characterized by the density matrix 

which represents a "phase-sensitive" state due to the presence of nonzero off-diagonal elements, thus al- 
lowing us to measure the phase properties of our setup. In fact, if the setup (cables and electronics) were 
ideal, the phase Z (glpg |e) — — Z (ej/Sq |g) associated with the off-diagonal elements of the matrix 



Pq of Eq. S7 would be zero. We can thus assume that any deviation from a zero phase corresponds to 
a phase error; third, we performed a single-qubit quantum state tomography (QST) on Qi, making possi- 
ble to measure experimentally /5q . Using our typical settings for a single-qubit QST', the time needed 
for each QST was approximately 8 s; fourth, we repeated a QST measurement every 8 s for a total time of 
40 minutes, corresponding to 300 measured density matrices; finally, we plotted Z (glpg^ |e) as a function of 



time. The so-obtained time trace is shown in Fig. SI \. Besides negligible slow-varying oscillations in the 



time trace [independent tests have shown that these oscillations might be due to temperature changes in the 
room-temperature cables (data not shown)], the overall histogram associated with the trace is approximately 
normally distributed about a mean value of Orad, with standard deviation tr^ ~ 0.065 rad (cf Fig. IS IB). 
However, we notice a general increase in the scatter of the time-trace data, as indicated by the dashed black 



lines in Fig. SI \. We thus divide the time trace in 10 sub-traces (time bins) with a time length of 4min 



each, compute the standard deviation for each sub-trace, and plot the so-obtained 10 standard deviations as 



a function of time. The result is displayed in Fig. SI 2, where the data is overlayed with a linear fit 



The plot of Fig. IS IC is useful in determining the phase errors associated with different types of two-qubit QST, 
as well as quantum process tomography (QPT) . In fact, two-qubit QST can be realized either by applying to 
each qubit the set of three unitary operations {/, Rx^^ , Ry^^} (I is the 2x2 identity matrix, Rx^^ a 7r/2 unitary 
rotation about the x-axis, and Ry^^ a 7r/2 unitary rotation about the y-axis), which we call "tomo," or the set of 
six unitary operations {/, Rx^^, Ry^^ ,Rx Ry Rx} (Rx '^^'^ is a — 7r/2 unitary rotation about the x-axis. 
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Ry ^^"^ a — 7r/2 unitary rotation about the y-axis, and i?J a tt unitary rotation about the .T-axis), which we call 
"octomo." 

In the case of two-qubit tomo, the number of operations that must be applied to the pair of qubits is given 
by the permutations of the allowed set of unitary operations, 3^ ~ 9. This number multiplied by the 4 possible 
joint probabilities for a two-qubit system, p^^, p^^, p^^, and p^^ (where, e.g., p^^ is the probability to measure the 
first qubit in the ground state with the second qubit in the excited state) gives a total of 36 probabilities. In the 
case of octomo, the total number of probabilities is given by the permutations of 6 unitary operations for 2 qubits, 
6^ = 36, times the 4 possible joint probabilities for a two-qubit system, for a total of 144 probabilities. 

In the experiments, the maximum likelihood estimator for each of the four probabilities p^„,p^^,p^„, and p^^ 
is obtained from the outcome of N measurements. We note that, in a joint two-qubit measurement each outcome 
consists of 4 numbers obtained simultaneously, where each number can be either or 1. The statistical sample 
consisting of N two-qubit joint measurements will be hereafter defined as Xj^, with l,m = g,e. Similarly 
to Eq. |S5] the maximum likelihood estimator (i.e., the estimated probability) for each of the four probabilities 
Pgg'Pge'Pee' ^^'^ Pee '^^us be obtained from 

1 ^ 

k=l 

where X^^^ represents the k-th outcome among the N measured. 

For a given fc, the four possible X^^, i.e., X^^, X^^, X^, and X^^, are measured simultaneously (with X^^ + 
+ Xj^^ + X^^ = 1). Hence, the effective number of events that has to be measured for each tomo is 36/4 = 9, 
and for each octomo 144/4 = 36. 

We typically measure 2500 events per second, and repeat each measurement N = 15000 times. As a conse- 
quence, a two-qubit tomo takes approximately 1 min, and a two-qubit octomo approximately 4 min. 

All data displayed in Fig. 2C of the main text were obtained using tomo, while all data in Fig. 3, D and E, were 
obtained using octomo. All density matrices used to reconstruct the x matrix of Fig. 3F in the main text were also 
obtained with octomo. The standard deviation due to phase errors can be estimated in each case by looking up the 
fit inFig.jST]::. 

Considering for example a two-qubit octomo with N = 15000, the statistical properties of the resulting density 
matrix p and of the corresponding metrics [fidelity J^, negativity J\f, concurrence C, and entanglement of formation 
£; cf. Ref.* and references therein for an extensive description of these metrics] are obtained as follows: 



(1) The probabilities P;„j associated with two-qubit octomo are estimated according to Eq. S8 As explained 
above, this corresponds to a total of 36 x 4 = 144 estimated probabilities. To simplify the notation, we will 
hereafter refer to these probabilities as P^, with i E {1,2,..., 144}; 

(2) The estimated probabilities P.^ are corrected for measurement errors [cf. Refs.^andl^for our standard pro- 
cedures to correct for measurement errors in the case of one and two qubits, respectively]. The corrected 
probabilities P^ are stored as a 144 x 1 column vector; 



(3) For each probability P^, the binomial standard deviation CTj, defined in Eq. S6 is calculated, thus obtaining, 
in the case of octomo, 144 different standard deviations; 

(4) For each of the 144 standard deviations CT[, calculated in (3), a set of M random numbers picked from a 
normal distribution with zero mean value and standard deviation (if, is generated. This results in a matrix of 
144 X M random numbers. Typically, M = 1000. 

By summing each column of such a matrix to the column vector containing the 144 estimated probabilities 
i-'j, we obtain a matrix of 144 x M probabilities, where each column simulates the result of a different QST 
experiment. 
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Fi gure S2: Confidence intervals for a density matrix and its metrics. (A) Histogram associated with the 6-th probability Pg of the 
vector of probabilities P^ , plotting the number of elements among the AI probabilities obtained in point (4) vs. the coiresponding value of the 
probability Pg . The data refers to the octomo for the state of Fig. 3E in the main text. The solid green line is a fit to a normal distribution. 
(B) Histogram for the real part of the matrix element with mean value (eg|p^ |ge) = 0.338 for the state of Fig. 3E in the main text. The 



For example. Fig. S2 \ shows the histogram associated with the 6-th probability Pg of the vector of proba- 



bilities Pj in the case of the octomo for the state of Fig. 3E in the main text; 

(5) Each column of the 144 x M matrix of probabilities obtained in point (4) is inverted by following the usual 
QST rules'^. This allows us to find the corresponding density matrix p™P''>'''^ vvith j e {1, 2, . . . , AI}, thus 
obtaining M density matrices associated with one state; 

(6) Physicality constraints are enforced on each, generally unphysical, density matrix p™^^'^^ by means of the 
MATLAB packages SeDuMi 1.21 and YALMIP (semi definite programming)'^. The physical constraints are 
such that each final - physical - density matrix should have unit trace and be positive semidefinite. 

In order to obtain the mean physical density matrix p associated with the M physical density matrices pj 
and the corresponding standard deviations, we calculate the mean value and standard deviation of the real 
and imaginary part of each matrix element for the M matrices pj. The mean physical matrix p will thus 
have elements {lm\p\pq) (with \pq) G A^2)' ^^ch of them (real and imaginary part) characterized by 
a given standard deviation. Figure |S2p shows the histogram for the real part of the matrix element with 
mean value (eg|/5^|ge) = 0.338 for the state p^ of Fig. 3E in the main text. As expected, the distribution is 
approximately Gaussian with a 0.95 confidence interval ±2a^^ = ± 0.005. 

The knowledge of the M matrices pj also allows us to estimate the confidence intervals for the relevant 
metrics characterizing the state p: T, TV, C, and E. This can easily be accomplished by calculating the 
metrics for each p^, thus obtaining M values for each metric, and then computing the mean value and 
standard deviation of the M values associated with each metric. 



We can follow a similar procedure to account for phase errors. We now pick two independent sets of M 
random numbers from a normal distribution with zero mean value and standard deviation (with opportunely 
estimated from Fig.|Sip depending on whether a tomo or octomo was used), thus generating two sets of M phase 

, Af}. In order to simulate phase errors acting independently on each qubit, 
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to a 4 X 4 measured density matrix p™'=''^, thus obtaining the j-th unphysical density matrix 

We can then proceed as in step (6) above and obtain a mean physical density matrix p and its statistical properties, 
as in the case of binomial-type errors. This allows us also to find the metrics associated with p and their statistical 



properties. Notice that the unitary transformation of Eq. S9 simulates random rotations along the z-axis of both 
qubit Q-^ and qubit Qj- 

The total mean physical density matrix is finally obtained by averaging the mean physical density matrix 
obtained in the case of binomial-type errors and the matrix obtained in the case of phase errors. The same applies 
to the mean values of all metrics. The corresponding standard deviations are found by summing in quadrature the 
values obtained in the case of binomial- type and phase errors. For example, the numerical value with confidence 
interval of each element of the density matrices in Fig. 2C of the main text were obtained following this procedure. 



These numbers are reported in Table S 1 



Incidentally, we found that phase errors do not add any significant contribution to the confidence intervals of 
the density matrix elements and of their metrics. 

Notice that, the reason why we decided to simulate the statistical properties of the data in Fig. 2 of the main 
text is because we only had 2 independent measurements of these data. Such a statistical ensemble is obviously 
insufficient to obtain reliable confidence intervals, which, thus, needed to be simulated. 



Experimental estimation of statistical errors 

In the case of the density matrices in Fig. 3E and of the Xm matrix of the quantum Fourier transform in Fig. 3F 
of the main text we had ensembles of independent measurements large enough to allow the confidence intervals 
estimation directly from the data. 

In particular, the density matrix /5q 28 in the left panel of Fig. 3E is the average of a statistical ensemble 
of M = 10 independent measurements, and the density matrices p^^^ and in the center and right panels 
of Fig. 3E, respectively, are the average of an ensemble of M = 6 independent measurements. The standard 
deviation of each matrix element (real and imaginary part) as well as the mean value and standard deviation of all 
metrics can easily be estimated from such statistical ensembles. 

Finally, the matrix Xm of Fig. 3F is the average of an ensemble of 15 independent measurements. This allows 
us to estimate the mean value and standard deviation of the process fidelity associated with the quantum 
Fourier transform (cf. main text). 



Statistical errors of fitted parameters 

The confidence intervals associated with the quantum phase tomography data shown in Fig. 4, C and F, of the 
main text are dominated by the statistical errors of the coefficients fitted from the data in Fig. 4, B and E, of the 
main text. In particular, the coefficient of interest is the phase of each curve in Fig. 4, B and E (or, more in general. 



of each curve in Fig. S12 C and D). 

We remind that the error vector associated with the vector of coefficients fitted to a curve is given by the square 
root of the vector S of the diagonal elements from the estimated covariance matrix of the coefficient estimates, 
{X"^ X)^^ (s)^. Here, X is the Jacobian of the fitted values with respect to the coefficients, X'^ is the transpose 
of X , and (s)^ is the mean squared error. This procedure allows us to estimate the errors associated with the 
fitted phases. These errors propagate through the quantum tomography process (cf. section on "Quantum phase 
tomography" at the end of these Methods), finally turning into the confidence intervals reported in Fig. 4, C and 
F, of the main text. 
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Definition of the qubit reference frame 

In this section, we briefly explain the concepts of reference frame and reference clock rate associated with a qubit. 
These concepts will be useful in understanding the dynamic phases acquired by the qubits when programming the 
quantum von Neumann architecture as well as the sequences used to tune up the CX-(j) gates and the XOR and M 
gate. 

In the two-level approximationE^l the Hamiltonian of a phase qubit can be written as 



/q(5) 

2 



HQ = h'-^a,, (Sll) 



with ground state |g) and excited state |e), and eigenenergies E„ and E^, respectively. In Eq. Sll /q(2) = 
AE^„/h ~ — £'g)//i represents the qubit transition frequency, which can be tuned by means of z-pulses with 
amplitude z, and is the usual spin 1 /2 Pauli operator. At the beginning of a CZ-(j) gate, each qubit is initialized 
in |g) at the so-called idle point, which corresponds to a z-pulse amplitude z = 0. The qubit transition frequency 
at the idle point is thus given by /q(z — 0) = Jq. 

In order to prepare a qubit in the excited state |e) or in a linear superposition |g) + |e), the qubit has to be 
driven by a microwave pulse. The Hamiltonian governing the interaction between the qubit and the microwave 
driving is given by ' ' 

i/o = hnj,{T)ay sin(27r/DT + , (S12) 

where ^^i''') '^he time-dependent driving amplitude expressed in unit hertz, /p the driving frequency, a-y the 
usual spin 1/2 Pauli operator, r the time, and ^^eiay arbitrary phase delay. By calibrating the microwave pulse 
such that the phase delay <j>^gi^y = 7r/2, we can rewrite the driving Hamiltonian as 

i7D = /ir!D(TU cos(2^/dt). (S13) 



By combining the qubit Hamiltonian of Ecj. |S 1 1 1 a nd the driving Hamiltonian of Eq. S13 we obtain the total 
Hamiltonian of the driven system, Hqj^ ~ Hq + H^. 

In our experiments the driving frequency /q is a fixed parameter that is set to be equal to the qubit transition 
frequency at the idle pointf^, 

/d = /q ■ 

For a given qubit, the microwave driving represents the reference frame associated with that qubit, with reference 
clock rate given by /q. Defining the detuning between the z-dependent qubit transition frequency /q(2;) and 

the reference clock rate /q as A(z) = /q(z) — /q, the qubit-driving Hamiltonian Hq^ can be expressed in the 

uniformly rotating reference frame by applying the unitary rotation D = e^^'^'^^^'^"'^'^^. The rotated Hamiltonian 
is thus given by 

Hq^ = D Hqo& -ihD-^D'f 

dr 

.-.M£),^+.^... (S14, 

where the counter-rotating terms have been already neglected. The dynamics associated with the pulse sequences 
used to tune up the CZ -cj) ga tes and the XOR and M gate can be understood by following the time-evolution of 

the Hamiltonian of Eq. S14 In particular, the Hamiltonian -ffgu describes the dynamic phases acquired by the 
qubits when they are brought outside their reference frame (i.e., qubit rotations about the z-axis). As it will appear 
clear when describing the tune-up sequences of the CZ-cf) gates and of the XOR and M gate, in the experiments 
we always compensate for such dynamic phases. 
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Programming the quantum von Neumann architecture 

The phase difference between the off-diagonal elements of the density matrices shown in Fig. 2C of the main text 
(red arrows) are due to the qubits being brought outside their reference frame during the pulse sequence in Fig. 2A 
(the qubits acquire dynamic phases), and to the angle accumulated by the microwave signal used to excite the 
qubits. The pulse sequence was calibrated such that the first density matrix has purely imaginary off-diagonal 
elements [cf. grey and overlayed red arrows in Fig. 2C (I) of the main text]. We can thus calculate the angles 
of the density matrices p^jj^ and p^y-^ by knowing the time duration of the various steps in the sequence and the 
corresponding qubit detunings (obtained from independent measurements), as shown by the grey arrows in the 
matrices of Fig. 2C, (III) and (V). As expected, the experimentally measured red arrows overlay the calculated 
grey arrows with high accuracy. We will later show a pulse method that allows us to compensate for dynamic 
phases during the experiment, rather than calibrating the phases a posteriori as in Fig. 2C. Such a compensation 
pulse method was used to implement the quantum Fourier transform and the XOR and M gate. 

The numerical values of all elements (real and imaginary part) of each density matrix in Fig. 2C of the main 
text are reported in Table S 1 . The confidence interval for the real and imaginary part of each complex number is 
also indicated. The confidence intervals correspond to two standard deviations (95 % confidence interval), where 
the standard deviations were calculated as explained in the section on "Statistical errors" of these Methods. 
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Figure S3: The CZ-(f> gate energy diagram. (A) Energy diagram for target qubit Qj^ coupled to bus resonator B, with Q^'s eigenstates 
indicated by |g). |e), and |f), and B's eigenstates indicated by |0) and |1). In general, states |fO) and |el) are detuned by a quantity <5g g. Thus, 

their effective Rabi frequency is given by (Sq g + 9q b)^'^^- Bloch sphere interpretation of the |fO)-|el) interaction. When <5g g = 0, 
Qj acquires a phase if> = tt (Top). When 5q g > 0, Qj acquires a phase < < tt (Bottom). 



The quantum Fourier transform 

The numerical values of all elements (real and imaginary part) of each density matrix in Fig. 3E of the main text 
are reported in Table S2. The 95 % confidence interval for the real and imaginary part of each complex number is 
also indicated. The standard deviations were calculated as explained in the section on "Statistical errors" of these 
Methods. 

As shown in the main text, the CZ-(/) gate is a fundamental element for the implementation of the quantum 
Fourier transform. In the rest of this section, we derive the analytical expression for the phase of a CZ-0 gate 
by diagonalizing the effective Hamiltonian of the Q]^-B-Q2 system and calculating its time evolution. We subse- 
quently describe the experimental pulse sequences required to tune up the CX-(f> gate and show three examples of 
Ramsey experiments used to measure the gate phase cf), when (f) = 0.01, = 7r/2, and (p — it. Finally, we discuss 
the origin of systematic errors in the measurement of the phase (j), showing that the global phase shift in the curve 
of Fig. 3C of the main text is due to a drift of the qubit operation point. 



Analytical expression of the phase of a CZ-c/j gate 

The CZ-0 gate demonstrated in the main text makes use of a bus resonator B that mediates the interaction between 
qubit Q]^ and Qj. During the CZ-0 gate qubit Qi is used as a qutrit, where the third eigenstate |f) plays an active 
role in the implementation of the gate. Qubit Qj^ represents the gate target and qubit Q2 the gate control. The 



energy diagram of the Qj^-B coupled system is displayed in Fig. S3 K. The coupled system consists of the states 
|g), |e), and |f) of the target qubit Qj^, and of the states |0) and |1) of the bus resonator B. In the rotating frame of 
resonator B and using the rotating-wave approximation, the system effective Hamiltonian can be written as 

= /iA|el)(el| + /i(A + 5^^^) |fO)(fO| + /i^(|el)(fO| + |fO)(el|) , (S15) 

where A represents the frequency detuning of with respect to the reference frame, 5q 3 the frequency detuning 
between states |fO) and |el), and g their on-resonance coupling. As already mentioned after Eq. S14 in the 
experiments we compensate the rotation about the z-axis of Qj^ associated with the detuning A. As a consequence, 



in Eq. S15 we can set A = and rewrite the system effective Hamiltonian as 



= hS^^j, |fO)(fO| + ;i^(|el)(f0| + |fO)(el|) . (S16) 
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The diagonalization of the Hamiltonian of Eq. S16 gives the eigenstates 



|QiB)_ = 

lQiB)+ = 
with eigenenergies 



1 + ^q,b/V'5q,b + .9q,b Iel> - ^1 - Sq,b/^%b + ~9q,b |fO> 

71 



(S17) 



(SI 8) 



E = 



QiB \/^Q^B+9q^B 



^QiB + y^Q^B + ffq^B 



(S19) 



(S20) 



respectively. 

Given the initial state |Q]^B)o = |el), after a time r the evolution U{t) = cxp{—iH^ffT/h) of the effective 
Hamiltonian of Eq. S16 acting on |Q]^B)q results in the state 



U{t) |QiB)o = ^ (\/ 1 + <5q^b/7^q,b + ~9iB lQiB>- e'^^- ^ 



+ \/1-'^q,b/V^q,b+.9q,bIQiB 



(S21) 



For a time corresponding to a full 27r-rotation between the states |fO) and |el), r — 1/{6q g + 5q b)^^^' Qi'B 
coupled system is in the state 



iQiB) = ^ [^1 + '^q,b/ V-^Ib + ~91b IQiB>- e-'('<5.^ -Z^"'^) 

+ ^1 - '5q,b/v/^C^^iqiB)+ e-^(^«.»^/^+'^); 

which, for simplicity, can be rewritten as 



where the phase (p is defined as 



iQiB) - e'^ |el) , 



'qiB + 9q^b 



(S22) 



(S23) 



(S24) 



Figure S3 3 depicts the Bloch sphere of the Q^-B coupled system for the states |el) and |fO), showing that the 
interaction dynamics between the two states always starts and ends at the same pole of the sphere. It is during this 
interaction that the phase (j) of Eq. S24 is acquired by the target qubit Qi. In particular, when (5q g = we obtain 
(j) = TT, while for i5q g > we obtain all phases < cf) < tt. 

An effective Hamiltonian similar to that of Eq. S16 governs the interaction dynamics for the Q2-B system. 
Hence, the interaction between the states of both the Qj^-B and Q2-B systems always starts and ends at the same 
pole of the coupling Bloch sphere. This has the important consequence that the CZ-ip gates used here are insensi- 
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Figure S4: Dynamic phase compensation for qubit Q-^. (A) Sequence without pulses on qubit Qj. Sequence steps: (la), Q-^ is initialized 
in state |Qi)jj^j = |g) at the idle point. The reference frame of Qj^ is indicated by a dash-dot magenta line. Resonator B, which is indicated 

by a dotted grey line, is in the vacuum state |0) ; (Ha), rotation .Rj}'^^ on Qj^ . The Gaussian pulse has a FWHM t^^^^ = 8 ns; (Ilia), z-pulse 
on Qj^ with amplitude Zj,^ ^ and length Tj,^ ^; (IVa), compensation pulse on Qj with amplitude z^-^^ and length t^^^\ (Va), rotation KjJ''^ 
on Qj^; (Via), measurement pulse on Qj^. (B) Sequence with pulses on qubit Qj. Sequence steps: (lb), the Qj^-Q2-B system is initialized 
in state |QiQ2B)jj(jj = |ggO), with both qubits at the idle point; (lib), rotation on Qj. The Gaussian pulse, in this case, has a FWHM 
Tp^jj^ = 7ns; (Illb), iSWAP between Qj and B with amplitude Zjg and length rjg = 24.97ns; (IVb)-(VIIlb), same as in steps (Ila)-(Vla) 
of A. Note that, in step (Vb) the z-pulse on Qj^ generates the phase <^> of the CZ-<ji gate. 

live to the relative phases of qubits Qj^ and Q2 when they are brought into resonance via resonator B. This feature 
allows us to use independent reference frames with incommensurate frequencies (and, hence, no special phase 
relationship) for each qubit, thus making possible to tune up each qubit with a separate calibration sequence. 
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CZ-0 gate tuneup 



Figure S4 A and B, shows the two sequences used to caHbrate the pulses applied to qubit during the CZ-(j> gate 
operation. We note that in the calibration sequence of Fig.[S4^ a series of pulses is applied to qubit Qi as well as 
to qubit Q2. We will show that by comparing the results obtained from the calibration of Qj^ without pulsing Q2 



(Fig. S4\ and Fig. S5 \) with those obtained by pulsing Q2 (Fig. |S4p and Fig. |S5p ) it is possible to measure the 
phase (p associated with the CZ-cf) gate (cf. also main text and Fig. 3C of the main text). 

Before delving into the analysis of the calibration sequences, we note that the idle point of qubits and Q2 
was set at a different position depending on the experiment. For the experiments of Fig. 2, C to E, in the main 
text, the idle point was set in between the memory and bus resonator for both qubits. This is also the case for the 
swap spectroscopies shown in Fig. IB of the main text. For all the other experiments, e.g., those described in this 
section, the idle point was set above the bus resonator for both qubits. 

Consistently with the vertical axis in Fig. IB of the main text, a z-pulse in the upward direction, which 
increases the qubit transition frequency, always corresponds to a negative z-pulse amplitude with respect to the 
qubit idle point. The opposite applies to the case of a z-pulse in the downward direction. 



The first calibration sequence for qubit Qj^, which is shown in Fig. S4 \, comprises the following steps 



(la) Qubit Qi is initialized in the ground state iQi)^;^) — |g) at the idle point, setting the qubit reference frame 



with reference clock rate /q . In Fig. S4 \, the reference frame is indicated by the dash-dot magenta line. 



During the entire calibration sequence, the bus resonator B is maintained in the vacuum state |0). Neverthe- 



less, in Fig. S4\ we indicate the presence of resonator B by a dotted grey line, which helps visualizing the 



frequency detuning between qubit and resonator; 

(Ila) Keeping the qubit detuning A = 0, a Gaussian microwave pulse with full width at half maximum (FWHM) 
■^FWHM is applied to Q^. The amplitude of th e pulse is chosen such that 27rri]-,T = 7r/2. In this case, the 
time evolution of the Hamiltonian of Eq. S14 yields a 7r/2 unitary rotation about the y-axisG3 Ry^"^, which 



brings the qubit into the new state |Qi)(ii^) = (|g) + |e))/v2; 

(Ilia) A z-pulse with amplitude z^.^^^, corresponding to a qubit frequency detuning A{zq2.-^) from the reference 
frame, brings the |e) |f) transition of Qi on or near resonance with B for a time Tp^-^- general, the 
z-pulse can be adjusted so that the |e) o |f) qubit transition is detuned by a frequency (5q g from resonator 



B (cf. Fig. S4 \, and also the section on the "Analytical expression for the phase of a CZ-ip gate" in these 
Methods, and the main text). During the z-pulse, since the |el)-|fO) transition of the Qj^-B coupled system is 
on or near resonance, the system remains always in the state iQj^B) = {\g) + \e))/y/2)^\0). This represents 
a dark state of the time evolution of the coupled system, yielding no swaps between and B, as opposed 
to the bright state iQjB) = (|g) + |e))/-\/2) Depending on the z-pulse amplitude ^cz-^ hence, 

on the detuning 5q g, the z-pulse length is chosen such that t^^-^ = 1/('^q b + .9q b)^^^ ('-f- section on 
the "Analytical expression for the phase (p of a CZ-(p gate" in these Methods). The z-pulse Zq.^ < moves Qi 



outside its reference frame. As a consequence, the time evolution of Eq. S14 acting on the state |Qi)(iia) 

gives the state |Qi)(i„,) = (|g) e-'"*-y» + |e) e+"^*")/%/2, where (b^y, = ~ A(^cz-0) This means 

that during the z-pulse Qj^ acquires an unwanted dynamic phase (j)^y„', 

(IVa) For the correct operation of the CZ-0 gate, the dynamic phase (j)^^^ must be compensated. This can be 
realized by applying a compensation z-pulse to Q^, with fixed length t^.^^ and variable amplitude z^mp- 
In order to avoid crossing the resonance with B, the amplit ude of the compensation pulse is swept in the 

opposite direction as compared to the z-pulse Zq^,.^ (cf. Fig. S4 \). In this case, the time evolution of i/go 
acting on |Qi)(m^) yields the state \Q^)^^^^^ = (|g) e"'*-/^ |e) e+^*./2)/y2, where 4/2 = - 

A(^cmp) ''cmp/2; 
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CZ-0.01 CZ-7I/2 CZ-n: 




-0.15 -0.10 -0.05 0.00 -0.15 -0.10 -0.05 0.00 -0.15 -0.10 -0.05 0.00 
cmp. amplitude, z (a.u.) cmp. amplitude, z (a.u.) cmp. amplitude, z (a.u.) 



Figure S5: Ramsey experiments for compensating the dynamic phase o f and measuring 0. (A) Probability of measuring Qj^ in 
|e), Pj, vs. compensation pulse amplitude z^^^ for the pulse sequence of Fig. |S4K. The compensation pulse is the z-pulse of step (IVa) in 
Fig.|S4j\. In the CZ-cf> gate experiments we always chose a compensation pulse length t^^^ = 7 ns. The blue dots represent measured data, 
while tne solid green lines are least-squares fits to a sine function. From left to right, the panels refer to a CZ-0.01 gate, a CZ-7r/2 gate, 
and a CZ-tt gate, respectively. In each panel, the vertical dotted black line indicates the amplitude z^.^^, chosen to compensate the dynamic 
phase The z^^^ numerical values expressed in the arbitrary units of our custom electronics are z^^^^ ~ — 0.135 for the CZ-0.01 gate, 
Zcmp — - 0.108 for the CZ-7r/2 gate, andz^^p ~ - 0.101 for the CZ-tt gate. Notice that the angle of the CZ-0.01 gate is different than that 
of the CZ-0.28 gate in the main text. This is because the CZ-0.01 gate shown here is the result of a single set of measurements, whereas the 
CZ-0.28 gate shown in the main text is the average of a set of 10 independent measurements. (B) Same as in A, but for the pulse sequence of 
Fig. |S4p . The relative phase between the Ramsey fringes of each panel in A and the corresponding panel in B gives the phase (j> of the CZ-ij> 
gate. For the three pairs of Ramsey fringes in this example, the relative phases are (j> = 0.01 rad, (f> = 7r/2 rad, and </> = Trrad. The vertical 
dotted black line in the rightmost panel is positioned at the same value of z^^^ as the corresponding panel in A, but, in this case, it indicates a 
minimum of the Ramsey fringe because of the tt shift introduced by the gate. 



(Va) A rotation Ry^^ similar to that in point (Ila) is applied to Qj^, bringing the qubit into the final state |Qi) (va) ~ 
-ism{cbj2)\g)+cos{cl)j2)\s)- 

(Via) Finally, a measiirement pulse is applied to in order to obtain the probability to find the qubit in |e), 
Pg = I cos{(j).^/2)\'^ = (1 + cos(j>^)/2. Since (p^ depends on the compensation pulse amplitude z^mp' '^he 
probability is also a function of z^.^^. In order to cancel the effect of the dynamic phase 0dy,j, ^cmp has to 
be chosen such that reaches a maximum, where the phase 0^ = 2Ktt, with K ^ Z. 
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In summary, the two qubit rotations Ry at the beginning and end of the calibration sequence of Fig. S4|\ 



reaHze a generalized Ramsey experiment, which allows us to measure the total phase acquired by Qj^ during the 



z-pulses that bring it outside its reference frame. The experimental data for the calibration sequence of Fig. S4\ 



are shown in Fig. S5 \, where the three Ramsey fringes are obtained for three different values of the detuning 



5q g, corresponding to a CZ-0.01, CZ-7r/2, and CZ-tt gate, respectively. For each Ramsey fringe in Fig. S5 \, the 



z-pulse amplitude z^^^ chosen to compensate the dynamic phase ip^^^ is indicated by a vertical dotted blacK line. 

The second calibration sequence for qubit Qj^ is shown in Fig.lS4b. The sequence, which is the same as the 
first calibration sequence with the addition of the pulses applied to qubit Q2, comprises the following steps: 

(lb) Qubit Q2, qubit Q2, and resonator B are initialized in the ground/vacuum state |Qi)(i[,) ® lQ2)(ib) ® l^)(ib) ~ 
|g) ® |g) ® |0), with both qubits biased at the idle point; 

(lib) A Gaussian microwave pulse with FWHM Tp^^^jy, is applied to Q2. The amplitude of the pulse is chosen 
such that fipT = TT. In this case, the time evolution of Hq^ acting on IQ2) fj,^-^ realizes a full qubit population 
transfer, Ry, bringing the qubit into the new state |Q2)(iib) = k); 

(Illb) The state |Q2)(nb) moved from Q2 to B by means of an iSWAP of length Tj^. At the end of the iSWAP, 
resonator B is in the state |B)(.jjjj,^ = 

(IVb) Qubit Q]^ is prepared in the state |Qi)(ivb) = (|g) + |e))/\/2 by means of a rotation i?^^^; 

(Vb) The same z-pulse as in point (Ilia) is applied to Q^. In this case, the Qj-B coupled system is in the state 
iQj^B) — (|g) + |e))/-\/2) ® which represents a bright state of the time evolution of the system, as 
opposed to the dark state IQ^B) = (|g) + |e))/\/2) ® |0). Depending 

'-'^ ^cz-</) ^'^'^ ^nd, thus, on ^q^B' 

at the end of the z-pulse the excited state |e) of Qj^ has acquired a phase [cf. main text; when 4> = n, the 
z-pulse is a — SWAP^ ^]. As in point (Ilia), at the end of such a pulse Qj^ has also acquired a dynamic phase 
resulting in the state iQi)^^,) = (|g) e^'-^'V, + |e) e+^^^^n e+»'^)/\/2; 

(VIb) The dynamic phase (f)^^^ acquired by Qi is compensated by means of a z-pulse, as in point (IVa). At the end 
of the compensation pulse Q^ is in the state |Qi)(vib) = [|g) e"'"*^'/^ + |e) e+'('*''/^+'^)]/\/2; 

(Vllb) A rotation Ry^"^ is applied to Q;^, bringing the qubit to the final state |Qi)(viib) ~ ~ ie"^/^ siii((/)a/2 + 
<^/2)|g)+e^*/2cos(0,/2 + 0/2)|e); 

(Vlllb) Finally, a measurement pulse is applied to Qi in order to obtain the probability to find the qubit in |e), 
= |e*'>/2 cos(0,/2 + 0/2) |2 = (1 + cos0b)/2, where cj)^ = <i>., + (f). 

The phase difference between the probability for the first and second calibration sequence allows us to 



measure the CZ-0 gate phase, (p^^ — (p^ ~ (p. This is illustrated in Fig. S5 where the phase difference between the 



Ramsey fringe in each panel of Fig. S5 \ and the corresponding fringe in each panel of Fig. |S5p gives the phase 
of a CZ-0.01, CZ-7r/2, and CZ-tt gate, respectively. 

The second calibration sequence can also be used to cross check the amplitude Z(,,^p of the compensation pulse 
chosen to cancel the dynamic phase (f)^^^- For example, when (p — n, the Ramsey fringe obtained from the second 
calibration sequence should reach a minimum for the same value of z^^^ for which it reaches a maximum in the 
first calibration sequence. This is confirmed by comparing the experimental data shown in the rightmost panel of 



Fig. S5\andFig.p5 



Figure S6 \ shows the sequence used to calibrate the pulses applied to qubit Q2 during the CTL-cp gate operation. 



The sequence comprises the following steps: 

(Ic) The system is initialized in the state |Qi)(i^) ® lQ2)(ic) ® l-^V^) ~ l^** ® Is) ® |0), with both qubits biased 
at the idle point; 
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Figure S6: Dynamic phase compensation for qubit Qj. (A) Sequence steps: (Ic), the Qj^-Qj-B system is initialized in state 
IQj^QjB)^!^^ = |ggO), with both qubits at the idle point. The reference frame of is indicated by a dash-dot magenta line. Resonator 

B is indicated by a dotted grey line; (lie), rotation Ry^^ on Qj. The Gaussian pulse has a FWHM Tp^Yfjjyj = 7ns; (IIIc), iSWAP between Qj 
and B with amplitude z^^ and length Tjj = 24.97ns; (IVc), z-pulse on Qj^ with amplitude ^^-^ length Tj,^.^; (Vc), iSWAP between B 

and Q2 with amplitude z-^^ and length Tjj = 24.97 ns. At the same time, compensation pulse on Qj^ with amplitude z^^^^ and length r^^p set 

in the sequence of Fig. S4B; (Vic), compensation pulse on with amplitude z^.„jp and length T;.,^p; (VIIc), rotation Ry^'^ on Qj; (VIIIc), 
measurement pulse on (B) Probability of measuring in |e), P^, vs. compensation pulse amplitude z^^^^ for the compensation pulse of 
step (Vic) in A. The blue dots represent measured data, while the solid green line is a least-squares fit to a sine function. The vertical dotted 
black line indicates the amplitude z^^^ chosen to compensate the dynamic phase acquired by Qj during the sequence in A, z^^^ ~ — 0.064. 

(lie) A rotation Ry^^ with FWHM Tp^pijy[ is applied to Q2; 
(IIIc) The state of Q2 is moved into B by means of an iSWAP; 

(IVc) Qi is moved on resonance or close to resonance with B through the same z-pulse as point (Vb); 

(Vc) The state of B is moved back to Q2 via an iSWAP. During and between the iS WAPs of point (IIIc) and (Vc), 
Q2 acquires an unwanted dynamic phase. 

At the same time as the iSWAP in point (Vc), the compensation pulse tuned up in the Qi calibration se- 



quence of Fig. S4 A and B [point (IVa) or (VIb)], is applied to Qi (this is not strictly necessary due to the 



independence of the calibration sequences for Qj^ and Q2); 
(Vic) A compensation pulse with fixed length t and variable amplitude z is applied to Q2; 
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Fi gure S7 : Qubit frequency drift. Phase (f) acquired by as a function of the detuning 6q q . The blue dots indicate the same experimental 
data as in Fig. 3C of the main text. The solid green line is the theory given by Eq. |S24| and the solid magenta line a fit to the function given by 
Eq.|S25| where 6^^,^^ is the only free fitting parameter. 



(VIIc) A rotation Rl^^ with FWHM Tp^j^,^ is applied to Q2; 



(VIIIc) The state |e) of Q2 is measured, thus obtaining the probabiHty as a function of z^^^ (cf. Fig. S6 3). 
Choosing a maximum of the probability allows us to cancel the effect of the unwanted dynamic phase 
acquired by Q2 during and between the two iSWAPs. 

Systematic errors 



The Ramsey fringes of Fig. S5 A and B, which are used to obtain the phase associated with three particular 
values of the detuning Sq g, can be extended to any arbitrary value of Sq g to obtain all possible values of (f). 
Figure 3C in the main text shows, for example, the phase cj) obtained for (5q g e [0, 70] MHz. In the figure, the 



theoretical expression for the phase (j) given by Eq. S24 (solid green line) is overlayed to the measured data (blue 



dots). The coupling gg g used to plot the theoretical curve was estimated from the time-domain swaps of Fig. 3B 
in the main text. A qualitative inspection of the figure shows that theory and measured data are shifted along the 
vertical axis by a detuning S^^.^f^, corresponding to an overall phase shift along the horizontal axis. 

The origin of the detuning S^^^f^, and of the corresponding phase shift, could be attributed to two main causes: 
(i) - Drift of the transition frequency of qubit Qj^ during the experiment; (ii) - drift in the room-temperature 
electronics. Each pair of Ramsey fringes used to obtain the phase cj) (e.g., the fringe in the leftmost panel of 



Fig. S5 \ and the corresponding fringe in the leftmost panel of Fig. [S5p) was measured within a few minutes. 
This time can be considered short enough to exclude the electronics drift as a main cause of the detuning between 
theory and data. We can thus assume the drift in the qubit transition frequency as the main reason for the detuning. 
This seems a fair assumption since the measurement of the time-domain swaps of Fig. 3B in the main text, which 
were used to calibrate the z-pulse amplitude Zcz-</) '^^'^ swapping time r^^ necessary to obtain each phase (j) 
(cf. main text and the two previous supporting sections on the CZ-cj) gate theory and tuneup), took approximately 
four hours. Both the swaps and Ramsey fringes were measured starting from large detuning (5q g = 70 MHz to 
zero detuning, resulting in a time delay between swaps and Ramsey fringes for each value of Sq g of approximately 
four hours. From independent measurements (not shown), in such a time interval we expect the qubit transition 
frequency to drift by a few mega hertz. 

In order to quantify the detuning 6^^^^^, we can fit the data of Fig. 3C in the main text with the function 

- ^^^^^ '^'i-ft (S25) 



The function of Eq. S25 represents the best fit we found for the data of Fig. 3C in the main text. Figure S7 shows 



the same data as Fig. 3C of the main text (blue dots), together with the theory of Eq. [S24| (solid green line), and 
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the fit of Eq. S25 (solid magenta line). The detuning obtained from the fit is \S^^^fi\ ~ 4 MHz, which is consistent 
with our expectation for a qubit frequency drift in approximately four hours. 

We note that the qubit frequency drift as well as the data scatter in Fig. 3C of the main text (or, equivalently, of 
Fig. S7 1 are peculiar to that measurement, where we intended to show all phases e (0, tt] in a single, long scan. 
When it will be required to use a specific phase to perform a quantum Fourier transform during an algorithm, 
we will first theoretically estimate the parameters z^^.^ and r^.^.^ and, then, search for the phase in the close 
vicinity of these parameters. This will allow us to measure only a small portion of the swaps of Fig. 3B in the 
main text and a few corresponding Ramsey fringes, which can be realized in a much shorter time than the long 
scan of Fig. 3C in the main text, thus significantly reducing the incidence of systematic errors. 



XOR gate and M gate tuneup 

In this section, we describe the experimental pulse sequences required to tune up the XOR gate and M gate shown 
in the main text. We also explain in detail the complete experimental sequence used to obtain one nontrivial 
entry of the truth table associated with the M gate. Finally, we outline the mathematical procedure at the basis of 
quantum phase tomography. 
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Figure S8: XOR gate tuneup. (A) (Left) Tuneup sequence (1-XOR). Qubit Qj is initialized in |g) at the idle point. Qubit Q2 and resonator 
B remain in state IQ2B) = |gO) during the whole sequence. The reference frame of Qj is indicated by a dash-dot magenta line. Resonator 
B is indicated by a dotted grey line. (Right) Ramsey fringe corresponding to the sequence on the left showing the probability of measuring 
Qj in |e), P^, plotted vs. z^^^^. The blue dots represent measured data, and the solid green line a least-squares fit to a sine function. The 
vertical dotted black line indicates the amplitude 2^.n,p, obtained from the fit, chosen to compensate the dynamic phase acquired by Qj^ during 
the first z-pulse in the sequence. (B) Tuneup sequence (2-XOR). As in A, but for Qj. (C) (Left) Tuneup sequence (3-XOR). The solid grey 
line indicates the ancilla qubit Qj, and the solid black line resonator B. The reference frame is set by Qj at the idle point, as indicated by the 
dash-dot magenta line. (Right) Ramsey fringe corresponding to the sequence on the left, where the probability of measuring in |e), P^, is 
plotted vs. Tj|.|. The blue dots represent measured data, and the solid green line a least-squares fit to a sine function. The vertical dotted black 
line indicates the time t^^^ , obtained from the fit, chosen to calibrate away the dynamic phase acquired by the state in B during and between 
the "write" and "read" iSWAR 
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XOR gate tuneup 



The quantum logic circuit of the XOR gate considered here is sketched in Fig. 4A of the main text. Figure S8 
A to C, shows the three sequences and the corresponding Ramsey fringes required to tune up the XOR gate. The 
concept behind each sequence is similar to the compensation of a dynamic phase acquired during a CZ-(j) gate, 
which has been elucidated in the section "The quantum Fourier transform" of these Methods. The third sequence 
needs special attention as it is applied to a resonator rather than a qubit state. 



(1-XOR) The first sequence, which is displayed in Fig. S8\ (Left), acts on control qubit Qj^. During the entire 
sequence, the control qubit Q2 and the target bus resonator B remain in the state IQ2B) ~ |gO), and all 
pulses acting on Q2 are turned off. Qubit is initialized in the ground state = |g) at the idle point. 

As for the case of the CZ-0 gate, the sequence consists of a Ramsey-type experiment used to determine and 
calibrate away the total dynamic phase acquired by during the XOR gate. Hence, the first step of the 
sequence is an Ry' unitary rotation on Q^, which brings the qubit to the equator of the Bloch sphere. 

Then, a z-pulse with amplitude %2-7r ^^'^ ''cz-tt — 39.08ns, brings the |e) -;-> |f) transition of Qi 
into resonance with B. The time Tq^.^^ is inversely proportional to the coupling strength between states 
iQj^B) — |el) and |fO). However, we note that in the sequence (1-XOR), B is always in the vacuum state 
|0) and, thus, no dynamics takes place between and B during the z-pulse. The only effect of the z-pulse 
is to detune outside its reference frame, which causes the qubit to acquire an unwanted dynamic phase 

</'dy,r 

In order to compensate (l>^y„, a second z-pulse must be applied to Qi. Such a pulse is characterized by 
a fixed time-length t^^^ = 5 ns and a variable amplitude z^^^. We note that all compensation pulses in 
the XOR- and M-gate tuneup sequences have the same time length of 5 ns. By continuously varying the 



amplitude Zcmc '^^e Ramsey fringe shown in Fig. |S8[\ (Right) is obtained. The dynamic phase (j)^^^^ is 



totally compensated when the probability reaches a maximum. In the figure, the vertical dotted black 
line indicates the compensation pulse amplitude chosen for this purpose, z^,^p ~ — 0.106; 



(2-XOR) The second tuneup sequence is displayed in Fig. |S8p (Left). The sequence is analogous to sequence (1- 
XOR), but acting on control qubit Q2 instead of In this case, z^^^ ~ - 0.127 [cf. Fig.lSSb (Right)]; 



(3-XOR) The third and last tuneup sequence of the XOR gate, which is shown in Fig. SB 2 (Left), acts on target 
resonator B. Sequence (3-XOR) represents a departure from the analogy between the tuneup of the XOR 
gate and the CZ-(j) gate, where only two compensation pulses were needed for the gate operation. 

As already explained in the main text, resonator B plays the role of the third qubit in our implementation 
of three-qubit phase gates. In order to use resonator B as an effective qubit, its state must be prepared and 
measured using either qubit or Q2 as an ancilla qubit. In our experiments, we have chosen Q2 to perform 
this function because of slightly better coherence times and measurement fidelities compared to Qj^. It is 
important to note that Q2 is actively used during the XOR gate. Consequently, the state of B can only be 
controlled before and/or after the gate operation. This issue does not constitute an experimental limitation 
since B represents the target of the gate and, thus, its state will not be controlled during the gate. However, 
once a state has been loaded in target B, it has to remain stored for a significantly longer time than any state 
stored in the control qubits and Qj. This is not an experimental limitation either, as the much longer 
coherence times of B compared to and Q2 (cf. caption of Fig. IB in the main text for numerical values) 
largely reduce the effect of a longer storing time. This experiment further proves the importance of the 
quantum von Neumann architecture, where the ability to store states in a memory makes possible to realize 
longer quantum computations. 

As for the control qubits Qj^ and Q2, also the state loaded in the target resonator B acquires a dynamic phase 
due to the detuning between the transition frequency of B and the reference clock rate of Q2. Note that only 
the frequency detuning with respect to Q2 contributes to the dynamic phase of the state in B because Q2 is 
the ancilla qubit chosen to manipulate and measure B. 
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Following the pulses in Fig. S8 Z (Left), the Ramsey expe riment necessary to compensate the dynamic 
phase acquired by B is indirectly preformed through Q^ r*^ ! *^ . Resonator B is initialized in the vacuum state 
|B) = |0), and qubit Qj^ and Q2 in the state IQ1Q2) = |gg), with both qubits at the idle point. While Qi 
remains in |g) during the whole sequence, Q2 is rotated into a linear superposition |g) + |e) by means of 
an Ry' rotation. Afterwards, an iSWAP with amplitude Zjg and time Tjg ~ 25.72ns is used to write the 
state from Q2 into B, which is thus prepared in the state |0) + This state remains loaded in B until it is 
read out by Q2 via a second iSWAP before the end of the sequence. Between the two iSWAPs, the control 
qubits Qi and Q2 remain in the state IQ1Q2) = |gg), and all pulses acting on both Qi and Q2 are turned 
off. Due to the excursion outside Q2's reference frame, the state |0) + |1) in B acquires a dynamic phase, 
which grows until the end of the readout iSWAP. Since the resonance frequency of B cannot be tuned, in 
order to calibrate away the effect of such a dynamic phase we delay the starting time of the readout iSWAP 
by a variable time t^^^. 

"'nil 

Finally, after the readout iSWAP, a second Ry rotation followed by a measurement pulse on Q2 completes 
the Ramsey experiment on B. The corresponding Ramsey fringe measured as a function of r^jg, is plotted 



in Fig. S8 Z (Right). Similar to sequence (1-XOR) and (2-XOR), choosing r^^] such that the Ramsey fringe 
reaches one maximum allows us to fully compensate the dynamic phase acquired by the state in B. The 
vertical dashed black line in the figure indicates the delay time chosen in the experiment, t^^^ ~ 3.46 ns. As 
a check, from the fit we also obtained a Ramsey fringe frequency of ~ 385.0 MHz, which agrees well with 
the Q2-B detuning ~ 369.8 MHz. 



M gate tuneup 



The quantum logic circuit of the M gate considered here is sketched in Fig. 4D of the main text. Figure S9 A to 
E, shows the five sequences and the corresponding Ramsey fringes required to tune up the M gate. The tuneup 
concept is similar to that used for the XOR gate, but with a few important differences due to the 1/2 CZ-tt gates. 



(1-M) The first tuneup sequence for the M gate, which is displayed in Fig. S9 \ (Left), acts on control qubit Q^. The 
only difference between this sequence and sequence (1-XOR) is that the single z-pulse with amplitude z^z,^ 
and length Tq2^_^ is now split into two z-pulses, z-pulse a and z-pulse /3, with amplitude cz-ir ~ -^cz-tt 
and time Ty^^^-Tr ~ '^CL-nl'^- This fact, however, does not affect the compensation of the dynamic phase 
acquired by Qj^, as the total excursion of Qj^ outside its reference frame remains unchanged. It is worth 
reminding that the z-pulse a and z-pulse /? bring the qubit |e) O |f) transition on resonance with bus 
resonator B. As a consequence, in sequence (1-M) the resonator remains in the vacuum state |0). 



From the Ramsey fringe in Fig. S9 \ (Right) we obtain one possible value of the compensation pulse am- 
pUtude that maximizes the probability of Q^, z^^^ ~ — 0.119; 



(2-M) The second tuneup sequence, which is displayed in Fig.|S9p (Left), acts on control qubit Q2. The sequence 
is the same as sequence (2-XOR) for the XOR gate. As in sequence (1-M), the z-pulse brings the qubit 
|e) o |f) transition on resonance with bus resonator B. 

From the Ramsey fringe in Fig. |S9p (Right) we obtain one possible value of the compensation pulse ampli- 
tude that maximizes the probability of Q2, z ~ — 0.077; 



(3-M) The third tuneup sequence, which is displayed in Fig. S9r (Left), acts again on control qubit Qj^, with 
control qubit Q2 and target resonator B either in state IQ2B) — |gO) or IQ2B) = |gl). In addition, all pulses 
acting on Q2 are turned off Qubit Qj^ is initialized in the ground state jQ^) — |g) at the idle point. 

In order to understand the dynamics of the interaction between Qj^ and B, we refer to the energy diagram of 



Fig.lSSK. If |B) = |0), after the first i?^^^ rotation on Q^, the Q^-B coupled system is in state IQ^) ® |B) = 
(|g) + |e))® |0). In this case, during the z-pulse a and z-pulse /? no dynamics takes place. As a consequence, 
in the time interval r^j^ between the two z-pulses, Qj^ remains biased at the idle point in state |Q]^) — |g) + |e). 
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Figure S9: M gate tuneup. (A) Sequence (1-M) for calibrating tlie i/zCZ-tt gates. (B) Same sequence as in Fig. [S8^. (C) (Left) 
Sequence (3-M) for compensating tlie dynamic pliase acquired by Qj during the shelving. Tj^.^: Time delay between the z-pulse a and 
z-pulse /3 (i.e., between the two i/2CZ-7r gates). (Right) Probability vs. rjj.|j^ for IQ2B) = |gO) (magenta squares) or IQ2B) = |gl) (blue 
dots). Dashed and solid green lines: Least-squares fit to data. (D) Same sequence as in Fig. |S8p . (E) As in C, but for bus resonator B. In this 
case, a detuning z^^^ is applied between the two 1/2 CZ-tt gates. No shelving is indicated by a solid magenta line, shelving by a hashed blue. 
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without acquiring any dynamic phase. The only dynamic phase acquired by Qj^ is that developed during 



the z-pulse a and z-pulse /3, which has already been compensated in sequence (1-M) (cf. Fig. S9\). The 
compensation pulse for such a dynamic phase remains turned on during sequence (3-M). 

If instead |B) = |1), after the first Ry^^ rotation on Q^, the Q^-B coupled system is in state \Qi) |B) = 
(|g) + |e)) (g) In this case, the reference clock rate is given by /q + /g. After the z-pulse a, the state 
iQj^B) — |el) gets shelved into the state IQ^B) = |fO) for a time t^j,, at the end of which the z-pulse (3 is 
applied. We remind that the |g) o |f) frequency is (2/q — (5,ji), where S^^i is the qubit nonlinearity defined 
as the frequency difference between the |e) o |f) and the |g) o |e) qubit transitions. In this experiment, 
/q^ ~ 7.2161GHz, (5„i ~ 140.6 MHz, and /g ~ 6.8150 GHz. During the time t^^, the coupled system is 
in state IQ^B) = (|gl) + |fO)) and acquires a dynamic phase (p^^^ — (/q — — /g) Tj,j,. This dynamic 
phase is independent from the dynamic phase acquired during the z-pulse a and z-pulse f3. 

In order to compensate the phase we delay the starting point of the z-pulse /3 by a time t^^.^- By 



continuously varying t^j^^, the two Ramsey fringes plotted in Fig. S9 Z (Right) are obtained. The magenta 
squares correspond to the case IQ2B) = |gO). As expected, in this case nothing happens as the dynamic 
phases acquired during the z-pulse a and z-pulse (3 were already corrected by the compensation pulse tuned 
up in sequence (1-M). The Ramsey fringe thus remains at the maximum probability chosen in that sequence. 
This fringe indicates that no resonator state has been shelved to the qutrit state \Qi) = |f). The blue dots, 
instead, correspond to the case IQ2B) = |gl). In this case, the sinusoidal dependence of the fringe clearly 
shows an excursion outside the Q^-B reference frame, which causes the dynamic phase (j)^^^ to be acquired 
by the state in Q^. This phase is totally compensated when the probability of Qi reaches a minimum. 
The reason why a minimum has to be chosen is because the Ramsey fringe is obtained with one excitation 
in the system, instead of no excitation as in sequences (1-M) and (2-M). This is analogous to the CZ-0 gate 
tuneup sequence displayed in Fig. |S4p . This can also be understood from the phase-gate cube of Fig. Slip 



(or the second row in Table S3 1, as the Ramsey fringe in sequence (3-M) measures the phase difference 
between vertex (5) and vertex (1) of the cube, which is Trrad instead of Orad. The vertical dotted black line 
in Fig.|S9p (Right) indicates the delay time chosen in the experiment, t^^^^ ~ 4.2 ns. As a check, a least- 
squares fit to the data (solid green line) allows us to extract the frequency of the Ramsey fringe, which is 
~ 248.2 MHz. This number is close to the expected value (/q — ^ni ~ /b) — 260.5 MHz (note that the data 
consists of one Ramsey oscillation period only. Hence, the ~ 10 MHz difference between the theoretically 
expected value and that obtained from the fit is within the fit confidence interval); 



(4-M) The fourth tuneup sequence, which is displayed in Fig.|S9p (Left), acts on target resonator B. The sequence 
is the same as sequence (3-XOR) for the XOR gate. In this sequence, the two iSWAPs bring the qubit 
|g) o |e) transition on resonance with bus resonator B. 

From the Ramsey fringe in Fig. |S9p (Right) we obtain one possible value of the delay time that maximizes 
the probability of Qj, t^^i2 — 2.77 ns (vertical dotted black line). Since the XOR gate and M gate 
experiments were performed several hours apart, the delay time for the M gate differs slightly from that for 
the XOR gate because of drifts in the qubit transition frequency (cf section on "Systematic errors" in these 
Supporting Online Material). In addition, from a least-squares fit to the data (solid green line) we extracted 
a Ramsey fringe frequency ~ 365.8 MHz, which agrees well with the detuning between the reference clock 
rate of Q2 and the transition frequency of B, ~ 369.8 MHz; 



(5-M) The fifth tuneup sequence, which is displayed in Fig.|S9p (Left), acts again on target resonator B. 

Sequence (3-M) served to compensate the dynamic phase acquired by during the shelving dynamics. 
Because B takes also part in the shelving, its state acquires a similar dynamic phase that must be compen- 
sated. Sequence (5-M) is analogous to sequence (3-M), with two differences. First, the Ramsey experiment 
is now performed on B via Q2 [as for sequence (4-M)]. Second, instead of adding a time delay, qubit Qi 
is detuned in the z direction by a quantity z^^^ for the entire interval between the z-pulse a and z-pulse 
(3. In fact, we are not allowed to use two times the same degree of freedom, i.e., the delay time Tj^u^ of 
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Figure S 10: M gate pulse sequence. Pulse sequence used to measure the Ramsey fringe in Fig. 4E of tlie main text (magenta dots). Tlie 
vertical dotted grey lines separate the sequence in 14 time frames. Each frame is described in the text. 



sequence (3-M), for the compensation of two independent dynamic phases. By continuously varying z^^^, 
the two Ramsey fringes plotted in Fig. |S9p (Right) are obtained. The interpretation of the fringes is the 
same as for sequence (3-M). The detuning value chosen in the experiment to compensate the dynamic phase 
acquired by the state in B during the shelving is z^^^ ~ — 0.012 (vertical dotted black line); 



(6-M) The sixth and last tuneup sequence for the M gate, which is not shown in Fig. S9 consists in repeating 
sequence (1-M). The compensation pulse for needs to be recalibrated due to the detuning z^g, set in 
sequence (5-M). The final value of the compensation pulse amplitude for Qj^ chosen in the experiment is 
^cmp^ -0.125. 

M gate pulse sequence 



Figure |STO| shows the complete pulse sequence utilized to measure the entry of the M gate truth table associated 
with state IQ1Q2) ® |B) = |ee) (8) (|0) + |1)). Step (1): Both control qubits Qj^ and Q2 and target resonator B are 
initialized in the ground state, jQj^QjB) = |ggO). The qubits are biased at the idle point. Step (2): Q2 is prepared 

in state |g) + |e) by means of an FCi^^ rotation. Step (3): The state in qubit Q2 is written into B via an iSWAR 
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Until this step, qubit Q2 serves as ancilla qubit to load resonator B. The iSWAP effectively zeros Q2, which can 
be now used as a control qubit in the M gate. Step (4): Both control qubits Qi and Q2 are loaded in state |e) by 
means of an Ry rotation. Step (5): First i/2CZ-7r gate between and B. Step (6): CZ-tt gate between Q2 and 
B. In the same time frame, the delay and detune necessary to compensate the dynamic phase on Qj^ and B due 
to the shelving are applied. Step (7): Second 1/2 CZ-tt gate between and B. Step (8): Compensation pulse on 
Q2. Step (9): Compensation pulse on Qi. Step (10): Compensation delay for the dynamic phase on B. Step (11): 
Zeroing gate applied to Q2. This step is necessary to re-use Q2 as ancilla qubit for controlling B. The zeroing 
is performed through an iSWAP between Q2 and Z2. Step (12): The state of B is read out by the zeroed Q2 via 
an iSWAP. Steps (13) and (14): A second Ry rotation on Q2 followed by a measurement pulse completes the 
Ramsey experiment on B. The Ramsey fringe obtained from this sequence is plotted in Fig. 4E of the main text 
(magenta dots). 



Quantum phase tomography 

The most general unitary operation describing a three-qubit controUed-phase quantum gate can be written as 
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(S26) 



In the ideal case, the amplitude of each diagonal element of Eq. S26 is unity, while the phase 



which state \lmn) e A^g is considered (cf. main text). All off-diagonal elements are zero. 



depends on 



It is straightforward to show that only seven of the eight phases of Eq. S26 are physically independent. In fact, 
the first complex exponential e* "^leso) can be factored out from the equation, allowing us to write the matrix 
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which is equivalent to [/"^ up to a global phase 1 



(S27) 
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The phases associated with each diagonal element in Eq. 



S27 



can be grouped in a column vector t defined as 
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with dimensions 



1) 



In the main text we have shown that by performing Ramsey experiments on the control qubits Qi and Q2 and 
on the target resonator B, it is possible to obtain the quantum phase tomography of the three-qubit XOR phase 
gate and of the Toffoli-class OR phase gate (M gate). In each Ramsey experiment one of the control qubits (or 
the target resonator) has to be prepared in a |g) + |e) (or |0) + |1)) state, while the other control qubit and the 
target resonator (or the two control qubits) are prepared in all four possible combinations of ground and excited 
state. In the case of an M gate, for example, the twelve states for each Ramsey experiment are reported in the first 
three columns of Table [S3] The fourth column shows the ideal value of the phase difference associated with each 
Ramsey experiment A similar Table can easily be obtained for the XOR gate (not shown). 

We note that the states of the control qubits and Q2 and of the target resonator B displayed in the first three 



columns of Table S3 constitute a general set of states for quantum phase tomography and, thus, can be used to 
characterize any type of three-qubit controlled-phase quantum gate. The phase differences associated with these 



Table S3: M gate Ramsey table. The first three columns indicate the state of the control qubits Qj and Qj and the state of the target 
resonator B for the twelve Ramsey experiments needed for quantum phase tomography. The fourth column shows the phase difference 
obtained from each Ramsey measurement for an ideal M gate. 
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states can be grouped in a column vector <p defined as 
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(S29) 



witii dimensions (12, 1). 

The aim of quantum phase tomography is to obtain the seven phase differences in vector r from the twelve 
phase differences in vector cp, which are measured by means of Ramsey experiments. In order to facilitate the 
explanation of quantum phase tomography, we now introduce a geometric representation of the phases associated 
with any three-qubit controlled-phase quantum gate. Figure |S1 1[ \ shows a cube, hereafter termed the phase- 
gate cube, the vertices of which contain information on the diagonal elements of a three-qubit controlled-phase 
quantum gate. The vertices of the phase-gate cube are enumerated according to the notation given in the space 
between Fig. 4C and Fig. 4F of the main text. Note that the phase-gate cube can directly be generalized to A^-qubit 
gates, in which case it has to be promoted to an A^-dimensional hypercube. 

Here, we will only consider the case of three-qubit gates, namely the XOR gate and M gate. The quantum 
logic circuits of these gates are shown in F ig. 4A and Fig. 4D of the main text, respectively. For convenience, 

[ B and D. The sign of each element t^°^, with k € {0, 1, . . . , 7}, of 
of vector r'^ (cf. main text) is given on the vertices of the phase-gate cube 



these circuits are also displayed in Fig 
vector T 



.XOR ^jjj each element rjf 



(cf. Fig. SI 1 C and E, respectively). As explained in the main text, a positive sign corresponds to a phase and a 



negative sign to a tt phase. As shown in Fig. S 1 1 C and E, the difference between the phases associated with each 



pair of vertices connected by a segment of the cube is indicated on the segment connecting that pair of vertices. 



For each gate, this gives a total of twelve phase differences corresponding to the elements in vector <p of Eq. S29 
We now use the phase-gate cube to determine the transformation matrix ^ between vector (p and vector 
T. Each row of T^^ must correspond to a segment of the phase-gate cube, and each column to a vertex, giving 
a matrix with dimensions (12, 8). The first row of T is associated with the phase difference between states 



|egO) and |ggO), 0|ego) ~ '^|ggO) ^'i- ^■^^ ^'^^ '■^^ ^^^^ °f '■^^ ^ g^te> Table S3 1. Adopting the enumeration 
in Fig. |S1 1| \, state |egO) corresponds to the vertex (4) of the phase-gate cube, and state |ggO) to the vertex (0). 
In the case of the M gate, for example, the first raw of T^^ must then be —1, 0, 0, 0, 1, 0, 0, 0. Following a 
similar procedure for all twelve segments of the phase-gate cube for the M gate, we readily find the entire M gate 
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transformation matrix T 
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(S30) 



Notably, the rank of the matrix T^^ of Eq. S30 is 7, as expected from the number of physically independent 
phases of the unitary matrix of a general three-qubit controlled-phase quantum gate, U"^ . A similar procedure can 
be used to obtain the transformation matrix associated with the XOR gate (not shown) or any other three-qubit 
controlled-phase quantum gate. 

Given the shape of vectors t and tp, and of the matrix ^, vectors t and are related by the simple linear 
system 



ip r 



T = if . 



(S31) 



Since the phase differences in cp are the only phases measured in the experiments, the system of Eq. |S31| has 
to be solved in order to find 



The matrix T is actually not invertible. However, the system of Eq. 



S31 



IS 

overconstrained by the experimental data and so it can be solved in a least-squares best fit sense, allowing us to 
obtain t . 

Figure |S12| \ shows the twelve Ramsey fringes used to measure the phase differences plotted in Fig. |S12p 
in the case of the XOR gate. Figure S12 B and D, shows similar results for the M gate. The phase differences 
associated with each Ramsey fringe are indicated in the space between panels A and B and, together with the 
corresponding pair of vertices of the phase-gate cube, in the space between panels C and D. Solving the system 



of Eq. S31 for the phase differences shown in Fig. S12 C and D, finally allows us to obtain the phases shown in 
Fig. 4, C and F, of the main text, thus realizing a full quantum phase tomography of the XOR and M gate. 

We note that the quantum phase tomography used here is inherently different from that developed in Ref.EH, 
where the time evolution of the quantum phase of the qubit state was used to infer information on the qubit 
dephasing mechanisms. 
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Figure S 1 1 : Geometric representation of the phases associated with the three-qubit XOR phase gate and the three-qubit ToffoU-cIass 
M gate. (A) The phase-gate cube. The eight vertices of the cube correspond to the diagonal elements of the gate unitary matrix. The vertices 
are numbered from (0) to (7), following the same enumeration as for the results of quantum phase tomography (cf. space between Fig. 4C and 
Fig. 4F in the main text). (B) Quantum logic circuit for the XOR gate, as in Fig. 4A of the main text. (C) The sign of each element t^'^^ of 
vector T™'' (cf. main text) is given on the vertices of the phase-gate cube. Each sign, +1 or —1, corresponds to a phase, or n (cf. main 
text). The difference between the phases associated with each pair of vertices connected by a segment of the cube (a total of twelve phase 
differences) is indicated on the segment connecting that pair of vertices. (D) Quantum logic circuit for the M gate, as in Fig. 4D of the main 
text. (E) As in C, but for the M gate. 
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Figure S12: Quantum phase tomography for the XOR and M gate. (A) Probability to measure either of the control qubits or 
Qj or target B in state |e) or |1) vs. Ramsey phase ip for the XOR gate. (B) As in A, but for the M gate. Open circles: Data. Solid lines: 
Least-squares fits to the data. The legend to the phases is indicated in the space between the panels. (C) The twelve phase differences in 
vector (p obtained from the Ramsey fringes in A. (D) As in C, but for the phase differences obtained from the Ramsey fringes in B. The pair 
of veilices of the phase-gate cube and the corresponding phase differences are indicated in the space between the panels. The error bars are 
due to the confidence intervals to the fits in A and B. Such confidence intervals propagate through the quantum phase tomography process 
generating the error bars in Fig. 4, C and F, of the main text. 
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